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Two bound, entangled fermions form a composite boson, which can be treated as an elementary 
boson as long as the Pauli principle does not affect the behavior of many such composite bosons. 
The departure of ideal bosonic behavior is quantified by the normalization ratio of multi-composite- 
boson states. We derive the two-fermion-states that extremize the normalization ratio for a fixed 
single-fermion purity P, and establish general tight bounds for this indicator. For very small purities, 
P < 1/N 2 , the upper and lower bounds converge, which allows to quantify accurately the departure 
from perfectly bosonic behavior, for any state of many composite bosons. 
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I. INTRODUCTION 

The composition principle, the treatment of composite 
particles as elementary objects despite their underlying 
structure, is a fundamental pillar of natural science pQ . In 
the microscopic world governed by quantum mechanics, 
from hadrons at the highest achievable energies [2 J to ul- 
tracold molecules [3] , the composition principle allows us 
to treat particles with integer spin as elementary bosons. 
This treatment can greatly simplify the understanding 
of the many-body behavior and the statistical physics of 
composite bosons; but when can we confidently apply it 
and reliably treat two bound fermions as a boson? 

The hierarchy of energy scales in nature would appear 
to provide an answer, and one intuitively expects that 
the bosonic behavior of bound fermions relies on their 
strong binding. For example, weakening the bound be- 
tween fermions indeed leads to the BEC-BCS crossover 
[3]. However, even when the constituents of cobosons 
[5], i.e. of compounds constituted of two fermions, are 
perfectly bound, it is not granted that the creation and 
annihilation operators of cobosons obey the bosonic com- 
mutation relations and exhibit perfect bosonic behavior: 
The Pauli principle for the underlying constituents may 
become relevant and thus jeopardize bosonic dynamics 
[5H71 I5HT5]. For good bosonic behavior, the occupation 
probability of any single-fermion state must be low, such 
that the constituent fermions of the cobosons do not com- 
pete for available single-fermion states. 

A satisfactory answer to the above question was given 
using the tools of quantum information [6j [7] : Indepen- 
dently of the actual physical system and of the binding 
strength between the constituents, the bosonic behav- 
ior of cobosons is intimately related to the entanglement 
between the constituting fermions [6j, and the impact 
of the Pauli principle fades away with increasing entan- 
glement [BJ IT3Hl6| . As an indicator for the entan- 
glement between the fermions, one may use the purity 
P of the reduced states of either fermion [17]: 1/P is 
the effective number of Schmidt modes, i.e. of populated 



single-fermion states. To treat cobosons as ideal bosons 
and evade the Pauli principle, there need to be many 
more single-particle states (1/P) than composites (A), 
i.e. N -P -C 1 [BJ. The original argument in [BJ was based 
on specific wavefunctions, it was generalized to arbitrary 
states in [7] , where general upper and lower bounds to the 
indicator of bosonic behavior - the ratio of normalization 
constants of many-coboson states that will be introduced 
below - were found. 

Here, we strengthen further the relationship between 
entanglement, as characterized by the purity of the 
single-fermion states, and the compositeness character 
of cobosons: We improve on existing bounds and derive 
the explicit form of those quantum states that maximize 
and minimize the normalization ratio for a given purity 
P. Our bounds are optimal, since they are saturated by 
the extremal states found. The tight upper bound comes 
close to the lower bound when A 2 • P <C 1, i.e. in this 
regime, not only is the deviation from perfectly bosonic 
behavior small, but it can also be bound very tightly via 
the purity. 

We first present the formalism for the treatment of 
A-coboson states and review previous results on the nor- 
malization ratio, in Section [IlJ following the notation of 
[BJ [7J. Our main result, tight bounds for the normal- 
ization ratio, is derived in Section |III| Examples and a 
discussion of limiting extremal cases are then given in 
Section [Tvl We conclude with a combinatorial interpreta- 
tion of the findings and an outlook on possible extensions 
and applications in Section [V] 

II. MANY-COBOSON STATES 

A. Normalization of many-coboson states and 
entanglement 

We follow the formalism of [BJ [7] for cobosons that 
are constituted of two fermions. The creation operator 
for a coboson constituted of distinguishable fermions can 
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always be written in the Schmidt decomposition [51 [7] 
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a, aJSJ, 



(i) 



i.e. as a sum over only one index, where the A-,- are the 

Schmidt coefficients, a] (Sj) creates an a (b) -type fermion 
in the Schmidt mode j, and the number of Schmidt coef- 
ficients is denoted by S. For cobosons composed of two 
indistinguishable fcrmions, we use the Slater decomposi- 
tion [8] 
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and set aj := flj-i ancl bj := /vL to recover the form Jl 



A state of N composite bosons reads 
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where the coboson normalization factor xn witnesses the 
possible departure from the familiar bosonic behavior 
[3 [H [111 EH- The factor X n, leading to the nor- 
malization of the many-coboson state, 1 = (N\N), is a 
completely symmetric polynomial in the Schmidt coeffi- 



cients Xj [20] 
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For ideal bosons, xn = 1 for all N, while = when 
the number of cobosons N is larger than the number of 
available fermionic single-particle states, S. 

The probability distribution A is characterized by its 
power-sums [151 |Tg| 120] 
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where M(l) = 1 due to normalization, < M(2) = P < 
1 is the purity of the distribution A. Power-sums are 
also called frequency moments, they are directly related 
to the Renyi entropy of the distribution A, H^ 6nyi (X) — 
log (M(m)) /(l — m). The M(m) are independent, but 
Jensen's and Holder's inequalities [3T] apply: 



M(k - 1)5=3 < M(k) < M(k - 



(6) 



Using power-sums, the normalization constant can be ex- 
pressed recursively [15] 
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where we set xo — 1 for convenience. 



Combinatorially speaking, the quantity xn is the prob- 
ability that, given a set of N objects which are each ran- 
domly assigned a property j (with 1 < j < S) with 
probability Xj, all objects carry different properties. For 
example, for 5* = 365 and Xj = 1/365, we obtain the 
solution to the "birthday problem", i.e. the likelihood 
that all members of a group of N people have a different 
birthday. The power-sum M(m) is the probability that, 
selecting m objects that each carry a property distributed 
according to Aj, all m objects have the same property. 
Therefore, we have the simple relationship \2 = 1 — M{2), 
while the xn with N > 3 are functions of all M(m) with 
in < N, as given by 



B. Normalization ratio as a measure of bosonic 
behavior 

The normalization factor xn of an iV-coboson state 
reflects the probability to create a state of N cobosons by 
the AT-fold application of the coboson creation operator 
on the vacuum. The resulting state reads 



(ct) 



A' 



10) 




(8) 



i.e. it is a weighted superposition of all possibilities to dis- 
tribute the pairs of fermions among the pairs of Schmidt 
modes (for species a and b). Every pair of Schmidt modes 
aj&j |0) is - at most - occupied by one fermion pair, in 
close analogy to the birthday problem. 

The normalization ratio, Xn+i/xn , has emerged as a 
decisive indicator for the bosonic behavior of a state of N 
cobosons: For example, adding an additional coboson to 
an iV-coboson state, i.e. applying the coboson creation 
operator, leads to the state [7] 
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i.e. the state is sub-normalized: It is possible that the 
addition of the (_/V+l)st coboson is inhibited by the Pauli 
principle, which occurs with probability 1 — Xn+i/xn ■ 
Similarly, the departure of the expectation value of the 
commutator [c, & \, which is unity in the ideal case, reads 



(jv|te,at]ijv} = 2^±i 

XN 



(10) 



The evaluation of xn scales prohibitively in the number 
of particles N, even when using the recursive formula 
([7]) [23]. Approximations to the normalization factor in 
terms of easily accessible quantities, such as the purity 
P = M(2), are thus desirable. From Eq. and for small 
N ■ P, a series expansion can be derived [131 [TS] , 
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1- N ■ P + N 2 (M(3) - P 2 ) 

+0 (N 3 (M(4) + 2P 3 - 2P M (3))) . (11) 
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On the other hand, an upper and a lower bound to the 
normalization ratio were found [7], 



l-P-N < 



Xn+i 
Xn 



< 1 -P. 



(12) 



However, the upper bound 1 — P is independent of N 
and cannot be saturated, and the form of typical states 
that maximize the ratio is not known. Here, we derive 
tight bounds, find the quantum states that saturate these 
bounds, and give a physical interpretation for their opti- 
mality. 



peaked distribution A^ with the same purity P = 1/5 
and 5 < 16. 

The peaked and the uniform distributions have ex- 
tremal properties: For example, they saturate the bounds 




III. TIGHT BOUNDS ON THE 
NORMALIZATION RATIO 

A. Extremal entangled states 

We are interested in the possible values of the normal- 
ization ratio Xn+i/x.n for a given P. In order to find the 
extremal values of xn+i/xn, we maximize and minimize 
this quantity under the constraints M(l) = 1, M(2) = P. 

Given a finite P, the number 5 of non-vanishing Xj is 
bound from below by L, the smallest integer that is equal 
to or larger than 1/P: 



S> L:= 



(13) 



FIG. 1: Visualization of probability distributions A with a 
maximal number of Schmidt coefficients S < 16 and purity 
P — 1/5. The diameter of each filled circle represents the 
probability Xj, its area is proportional to A^. The diameters 
sum to unity, while the total area adds up to P. Distributions 
with the same occupied gray area thus have the same purity 
P, the three distributions shown cannot be distinguished via 
their purity P, but only through higher-order power-sums. 
Given P — 1/5, different distributions A can lead to different 
normalization factors \n- The uniform distribution X^ min- 
imizes the normalization ratio, while the peaked distribution 
A* p ' maximizes it under the chosen constraint S < 16. 



on the higher-order power-sums M(k > 3), given by (JsJ) . 
In the limit 5 — > oo, we have 



Distributions A with 5—1 equal coefficients [35] con- 
stitute extremal states, and will turn out to minimize or 
maximize Xn+i/xn- We thus define 

. (±) l± v /(5-l)(5P-l) 
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When we choose 5 = L, we obtain the uniform distribu- 
tion \M := X ( (s l L) with A^ } < X [ ^ {2 Jj} . This distri- 
bution minimizes the number of non-vanishing Schmidt 
coefficients. 

For 5 > L, we find a peaked distribution XS P ^ 

that satisfies A^ 



peaked coefficient converges to VP, while all other 
coefficients become vanishingly small, while the distribu- 
tion always remains normalized and possesses the purity 
P. Choosing P = 1/5 implies L — S and A^ = A^ - 
all coefficients are then identical, a maximally entangled 
state is obtained. In Fig. [I] we illustrate the uniform dis- 
tribution A( u ) , a randomly chosen distribution A, and the 
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(S>L) 

In the limit 5 — > oo, the 



(15) 



and, for fractional values P = l/£, we find for the uni- 
form distribution: 



M {u) {k) = P-M (u \k-l). 



(16) 



Combinatorially speaking, given the probability to find 
a pair of objects with the same property {i.e. given the 
value P = M(2)), the probability M(3) to find three 
objects with the same property for three randomly chosen 
objects is maximized by the X^ and minimized by the 
\( u ) distributions. 



B. Normalization ratio for extremal states 



Since only two different non- vanishing values of Xj ap- 
pear for the uniform and the peaked distributions, the 
normalization ratio can be computed explicitly for these 
extremal distributions: 
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X ( n ] (L - 1) (i + £(JV - 1)P - JV (l + y/(L - l)(LP - 1))) 



(17) 
(18) 



where Xjv is the normalization factor for the uniform/peaked distribution. We can now formulate our tight bounds 
on Xn+i/xn, given an arbitrary distribution A of 5 Schmidt coefficients: 

1-pV^4^ ( <^ ( <%< lim % ( ^^(P)<1-P, (19) 



where Xn+i/Xn 1S computed for the finite 5 defined by 
A, and we define the upper bound 

P-N 



U N (P) = 1 



1 + (JV- I) VP 



(20) 



The inequalities ( 12 ) are represented here by the ex- 



tremal lower and upper bounds (i,v), they were first 
shown in [7], for an alternative proof see [19]. We prove 
the new bounds (ii — iv) in Appendix [A] and discuss their 
physical implications in the following Section |TV} 



When the purity is decreased for a constant 5, the up- 
per and lower bounds eventually merge when P = 1/5 is 
attained (see also inset (a)). 

The upper bound can be expanded in powers of vP: 



u N (p) = i-J2 p 
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N((N - l)P + (1 - AT)«P(™+ 1 )/ 2 ) 



(jV-l)(l + (iV-l)VP) 



IV. ILLUSTRATION AND INTERPRETATION 
OF BOUNDS 

All bounds for N — 2 are illustrated in Fig. [2j 



A. Upper and lower bounds 

The authors of [7 showed that the lower bound 
1 — N ■ P in (12 1 is attained for fractional values of P, 
i.e. for P = ljL = 1/S. Setting L = 1/P in Q re- 
produces the bound, and Xjv+i /x { n ] = 1 - iV • P. The 



saturation can also be observed in Fig. [2] The tight lower 
bound (blue line) coincides with 1 — N ■ P (red line) for 
fractional values of the purity P = 1/S*. When P -C 1, 
and thus P « 1/L, the tight bound (18) differs only 
marginally, while for large purities P ^ 1/N the tight 
bound can differ significantly from 1 — N ■ P. In the limit 
P — » 1/iV, the tight bound (18) and the previous bound 
1 - N ■ P differ by a factor (1 + iV) /2. 

In contrast to the previously established upper bound 
1 — P, our tight upper bound ( |19| w) depends on the num- 
ber of particles N. When the number of non-vanishing 
Schmidt coefficients S is finite, the bound |l7] ) is more 
efficient than the limiting case S — > oo: In Fig. [2j the 
dashed black lines show the upper bound for finite S, 
while the blue line indicates the absolute upper bound. 



with convergence radius P < 1/(N — l) 2 . To second 
order, we then find 

U N {P) ml- P-N + P 3 / 2 (N 2 -N) + C(P 2 ), (22) 

i.e. the upper and lower bounds coincide in the limit 
P -)• 0. Indeed, for P < l/(N- l) 2 « 1/N 2 , the denom- 
inator in (20), 1 + (JV — 1)VP, is of the order of unity. 
This behavior is illustrated in Fig. [3| where we plot the 
deviation from perfect bosonic behavior, 1 — xn+i/xn- 
The JV-dependence of the new upper bound is apparent, 
as well as the convergence of the upper bound to the 
lower bound. In particular, the purity essentially defines 
the normalization ratio in the range P • JV 2 <C 1 . 

In order to illustrate the typical behavior of random co- 
bosons, we generated 3- 10 s random distributions A [2"7] . 
sampled according to the Haar-measure [27H29] , Pairs 
(P, Xs/Xi) are counted in a grid with 1000x1000 bins, 
which is translated to a color-code in Fig. [2] We gen- 
erated states with 5 = 3 (main figure), 5 = 4 (inset 
(b)) and 5 = 5 (inset (c)) non- vanishing Schmidt coef- 
ficients. The bounds for finite 5 are indeed reached by 
randomly generated states. We also observe a concentra- 
tion of states around the peak value of (P, X3/X2) when 
the number of Schmidt modes 5 is increased: The vast 
majority of randomly generated states in high dimensions 
share very similar entanglement properties [30] . 
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Counts per bin 




FIG. 2: (color online) Bounds for the normalization ratio 
X3 /X2 as a function of the purity P. The red solid lines indi- 
cate the extremal bounds in Eq. ( |19[ ) found in 0- The solid 
blue lines denote the tight bounds in Eq. ( 19 1 that can be 
achieved for any value of P. For a given maximal number of 
Schmidt-coefficients S, the black dashed lines are the corre- 
sponding upper bounds (for peaked states X^), which merge 
with the lower bound at P — 1/S . Inset (a): The scale in 
P is logarithmic, the (black dashed) upper bound for S = 50 
is close to the total upper bound at P = 1/10, but it merges 
with the lower bound at P = 1/50. The distributions show 
the numerically obtained values for the normalization factor 
and the purity, for different fixed Schmidt numbers, S = 3 
(main figure) and S — 4, 5 (insets (b), (c)). 



B. Limit of many particles 

Surprisingly, increasing the number of particles at con- 
stant purity P does not always fully destroy bosonic be- 
havior: The lower bound 1 — TV ■ P admittedly decreases 
with TV when P > and vanishes for P — 1/TV - the cor- 
responding uniform state X^ u ' consists of a finite number 
L = [~1/P] of Schmidt modes, such that at most L par- 
ticles can be accommodated and xl+i = 0. The peaked 
state, however, leads to non- vanishing xn+i/xn for ar- 
bitrarily large particle numbers: In the limit TV — > oo, we 
have 



U N {P) N A°° 1 - 



P. 



(23) 



That is to say, for the peaked state the departure from 
bosonic behavior, as quantified by the ratio \n+i/xn, 
amounts to at most VP, for any number of particles TV. 
This counter-intuitive result can be understood by the 
extremal form of the distribution X&': When a fermion 
pair populates a Schmidt mode other than the first one 



(which is populated with probability VP), it can essen- 
tially be neglected for the impact of the Pauli principle 
on the next fermion pair, since there are arbitrarily many 
such modes that are occupied with vanishing probability 
(in the limit S — > oo). Assuming that TV 3> 1 particles 
were successfully prepared, the probability that the first 
pair of Schmidt modes is populated by some fermion pair 
/ — N 

is 1 — VP , i.e. very close to unity. Adding an (TV + l)st 
particle is thus successful when this last particle does not 
end in the first Schmidt mod e, i. e. the success probability 
is 1 — VP, just as given in (23). Colloquially speaking, 



there are always enough Schmidt modes to accommodate 
another particle. The last added particle must not, how- 
ever, end in the first Schmidt mode, since the latter is 
occupied with nearly unit probability. 

Although an TV-coboson state still behaves bosonic to a 
certain degree, the normalization factor for such a state, 
given the peaked distribution , is 

X ^ = (i-Vp) N ~\i + (N-1)Vp), (24) 

which converges to zero for any < P < 1 in the limit 
TV — > oo. The analogous normalization factor for a uni- 
form distribution (assuming P = 1/L) reads 
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(25) 



which decays faster in N than Xn an< ^ vanishes identi- 
cally for N > L. 
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FIG. 3: (color online) Deviation of the normalization ratio 
from unity, 1 — Xn+i/xn, as a function of the purity P, in log- 
log-representation. The red dashed line is the TV-independent 
upper bound to xn+i/xn, i.e. here it represents the lower 
bound P to the deviation. The black dot-dashed, green solid 
and blue dotted lines indicate the tight bounds for N = 1000, 
N = 100, and TV = 10, respectively. The previously found 
bound 1 - TV • P and the tight bound (fel do not differ 
significantly in the present regime P < 1, they thus cannot 
be distinguished in the plot. The normalization ratio can take 
any value in the respective shaded areas. The deviation from 
ideal behavior as well as the gap between the upper and lower 
bound decrease with decreasing P. 
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V. CONCLUSIONS AND OUTLOOK 

The A^-coboson normalization ratio Xn+i/xn was es- 
tablished as an important indicator for the composite 
behavior of non-elementary bosons [BJ. Our main result 
is a new, A^-dependent upper bound for xn+i/xn, and 
the explicit representation of the states that saturate this 
bound. 

In the limit of small purities P < 1/N 2 , the bosonic 
behavior of an A^-coboson state is very tightly defined 
by P, since the lower and upper bounds merge. In prac- 
tice, the purity of a bound pair of particles often satisfies 
P < 1/Af 2 , e.g. for atoms in a trapped BEC pj ITU] . 
Our bounds thus provide a simple and reliable way to 
quickly check the departure of bosonic behavior of any 
type of cobosons. With a combinatorial argument, we 
can understand this clear determination of bosonic be- 
havior: For small purities PN 2 < 1, the probability to 
not finding A~ objects with different properties is essen- 
tially determined by the probability to find exactly one 
pair of objects with the same property, which is defined 
by M(2). Triplets and larger combinations that depend 
on higher-order M{k) are then essentially negligible. 

When the purity is not very small, P ~ 1/A~, however, 
the form of the wavefunction does play a role for bosonic 
behavior, which is then not entirely defined by P. It 
might be possible to access entanglement properties of 
bound fermions via the higher power-sums JVf(m), which 
can be obtained by measuring xn [31] . 

Formally speaking, we found bounds to the completely 
symmetric polynomial Q in terms of the first and sec- 
ond power-sums ((5|, M[l) = 1 and M(2) = P. For 
tighter bounds, one could specify also the third power- 
sum, M(3), and repeat the maximizing- and minimizing 
procedure of Section [X] to find those states that min- 
imize/maximize Xn+i/xn for given M(l), M(2), M(3); 
this procedure could be extended to even higher orders. 
It is, however, not immediate how operations that are 
analogous to (A3 1 and (A4) but which leave M(3) in- 
variant can be constructed. In the typically encountered 
domain of small purities, P <C this endeavor is not 

an urgent desideratum, since the encountered bounds as 
a function of P are already tight. On the other hand, 
using relations between Renyi entropies of different or- 
ders [221 [53] , our bounds can be re-formulated in terms 
of other indicators for entanglement, such as the Shannon 
entropy of A. 

Given a fixed purity P, the uniform distribution X^ 
minimizes the probability that the Pauli principle is ir- 
relevant, while the peaked distribution maximizes it. In 
other words, the A^-coefficient (M (3) — P 2 ) in the expan- 



sion (111 is maximized. Although a peaked or canyon dis- 



tribution leads, in general, to a normalization ratio that 
is smaller than for the uniform distribution [15], this is 
mainly due to the consequent change of purity. For fixed 
purity P, the bosonic behavior of the uniform distribu- 
tion is actually inferior with respect to the peaked one. 
Combinatorially speaking, we have considered a vari- 



ant of the birthday problem with non-uniform probabili- 
ties [24] . Here, Schmidt modes or single-particle quantum 
states take the role of birthdays [24, 2BJ or surnames |32j . 
Rather counter-intuitively, the optimal bosonic behavior 
for a fixed purity P = M(2) is found by maximizing the 
probability to find three objects with the same proper- 
ties, i.e. M(3). This result can be understood as follows: 
Any pair of objects that have the same property is as 
deleterious as any triplet (or any other m-tuplet). The 
probability to find a pair, however, decreases with in- 
creasing M(3). This decrease has a larger impact on the 
overall probability to find all objects with different prop- 
erties than the consequent increase of the probability to 
find triplets with M(3). For example, for A" = 3, the 
probability to find a pair amounts to 3(P — M(3)), the 
probability to find a triplet is M(3). The overall proba- 
bility to find all objects with different properties amounts 
to 1-3P + 2M(3). 

To complement the analytical bounds, we have nu- 
merically generated random states, which do not only 
show a concentration around the most probable value of 
the purity P [3Uj, but they also cluster around a cer- 
tain value of the normalization ratio, consistent with the 
concentration-of-measure phenomenon. It remains to be 
studied how random states in higher dimensions behave 
in general, i.e. what is the typical normalization ratio 
and the distribution of states with a given purity. 
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Appendix A: Proof of bounds and tightness 



In order to show Eq. (19) (ii) and (Hi), we co nstruct 

These 



Al 



operations on the distribution A, in Section 
leave the sum of the Xj and the purity P invariant, 
while they increase or decrease the normalization ratio, 
as shown in Section A 2 Since only the extremal distri- 
butions A^ and A^ remain invariant under the applica- 
tion of the operations, these distributions maximize and 
minimize the normalization ratio, respectively. 



1. Uniforming and peaking operations 

We construct uniforming and peaking operations on 
the distributions A that act only on three selected A,,-, 
with the indices 1 < ji < 32 < J3 < S. The operations 
will leave 



K\ — Xj 1 + Xj 2 + Xj 3 , 
K 2 = X 2 h +X 2 2 +X] 3 , 



(Al) 
(A2) 
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invariant, and, consequently, also the total sum of the Xj conjecture 
and of the A?. The third power-sum, M(3), however, will 
be changed by the operations. 

In analogy to an analysis of the birthday-problem with 
non- uniform birthday probabilities |26j , we define the two 
operations, T u and T p , on the probability distribution A: 



TP /u {Xh) = T P/u {Xh) = 



32,33) 



1 
3 
1 

6 



Ki ± ^J6K 2 - 2K\ 
2K X t J6K 2 - 2Kf 



(A3) 



where the upper (lower) sign in ± and =F refers to the 
peaking (uniforming) operation T p (T u ). For K\ < 2K 2 , 
we formally have F u (Aj 1 ) < 0, and we alternatively set 



r»(A ix ) = 0, 



(A4) 



For convenience of notation, we set 

= r"(Aj), X P = T p (Xj). (A5) 

Colloquially speaking, T u levels out the three coeffi- 
cients Xj 1 , Xj 2 , Xj 3 and thus makes the distribution more 
uniform, whereas T p makes the distribution more peaked. 
In both cases, the purity P is kept constant. The oper- 
ations push a distribution A towards the uniform and 
peaked distribution, respectively, as illustrated in Fig. [4] 




(A£,A£,A£) (A J15 X h , X h ) (A?, X p h , X P J3 ) 

FIG. 4: Action of uniforming and peaking operations T u and 
r p . The tuple (Xj 1 , Xj 2 , Xj 3 ) is leveled out by T u , and made 
more peaked by T p . We show only the coefficients Xj with 
indices ji,j2,j3, all other coefficients At, remain constant un- 
der the application of the operations, for the choice of indices 
iijj2,j3. The gray area represents the sum of the squared 
coefficients, it is the same for all three distributions. 

The uniform (peaked) distribution X^ u ' (A^) is the 
only one that remains invariant under the application of 
pp (F M ), for all choices of ji, (disregarding permu- 
tations of the indices), which can be seen by applying the 
operations on the distributions. 



2. Normalization ratio under operations 

We now show that the uniforming (peaking) operation 
p u (p) reduces (increases) the normalization ratio, i.e. we 



X iv+i(r"(A)) < xn+i{x) < X N+i(r p (X)) 

Xiv(r«(A)) " X n(X) ~ Xn(T p (X)) 



(A6) 



where we made the dependence of xn on the distributions 
explicit. 

Without restrictions of generality and for convenience 
of notation, wc set j\ = h J2 = 2,j3 = 3, i.e. we let the 
operations act on the first three Schmidt coefficients. 

We define 



Xn — Xn(Xa,---,Xs), 
A = A1A2A3 



(A7) 
(A8) 



i.e. formally xn is a normalization factor, but for an 
unnormalized distribution {A4, . . . , A5}. With these def- 
initions, we rewrite xn as 

Xn — A ■ xn-3 + (A1A2 + A3A2 + AiA 3 )xjv-2 

+(A 1 +A 2 + A 3 )xAr-i+XJV- (A9) 



The terms 



AiA 2 +A 3 A 2 + AiA 3 = I (K\ - K 2 ) , (AI0) 



Ai+A 2 + A 3 = K u 



(All) 



and Xk with k G {N — 3, ... , N} do not change upon 
application of r u / p , i.e. only the product A is affected by 
the operations. 

Conjecture: It holds 



A" A2A3 < A < X P X P X P , 



(A12) 



where X u ^ p is the result of the operation T u / p on Aj . 

Proof: We re- write the products in terms of K\ , K 2 
and Xi 



A l A 2 A 3 

A 

X\ X 2 X P 



L (k x - y/6K 2 - 2Kfj ioxKl>2K 2 
(2Ki + ^/6K 2 - 2Kfj 
for Kf < 2K 2 



108 

X 



-Ai (2X\-2X l K 1 +Kl-K 2 



(A13) 



108 



K, + J6K 2 - 2K{ 



x[2K 1 - J6K 2 - 2K{ 



For given K\ and K 2 , we can find the possible A2/3, 
leaving Ai as a free parameter, 



1 



A2/3 = -[K 1 -Xi± yJ2XiKi - Kf - 3X1 + 2K 2 
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Ai' 



Al,: 



\P \P \P 



A l A 2 A 3 




Ai 



Al 5 maxf 



FIG. 5: (color online) Upper panel: Possible values of Ai, A2 
and A3, given K\ and K2, as a function of Ai, in arbitrary 
units. Lower panel: Behavior of A1A2A3 as a function of Ai, 
and lower and upper bounds. The lower (upper) bounds are 
attained if and only if the configuration fulfills \j 1 > A J2 = 
Xj 3 (Aj! < A J2 = A J3 ), where {ji, 02,33} is a permutation 
of {f , 2, 3}. These points are marked by thin dotted vertical 
lines. 



The requirement A2/3 > gives 
1 



Ai. 



max/min 
•^l.min 



Kt ± V2J3K 2 -iq), (A14) 



< Ai < Ai 



(A15) 



Give n K \,Ki, all possible values of Ai then fulfill 
Eq. ((AT2l. □ 



The inequalities in (A12| are saturated for A 



:u/ P 



A, 



(modulo permutation of the indices). Possible values of 
Xj and the behavior of the product A1A2A3 = A are shown 
in Fig. [5] With ( A9 ) , it immediately follows that 



( u ) ^ ^ (p) 
Xn < Xn < Xn ■ 



(A16) 



Using ( A12 1, we can set 



AfA^Af =: A(l + e), 
A" A2A3 =:A(l-6), 



with e, S > 0, and 
Bn = 



(A1A2 + A3A2 + AiA3)xat_2 

+(Ai + A 2 + A 3 )xjv-i + Xn- 



For any distribution A (which does not need to fulfill 
J2jXj — 1), the Newton-Maclaurin inequality holds J7J 
20] . which reads 



Xn+i 



< 



XN 



Xn Xn-i 

We thus have 

Xn~3 > Xjv-2, Bn > Bn + i. 



(A17) 



(A18) 



Our original conjecture (A6) is equivalent to 



A(l - S)xn-2 + B N+ i Axn-2 + B N+ i 



A(l - S)xn- 3 + B 



N 



Axat-3 + B 



N 



< 



A(l + e)xN-2 + gjv+i 
A(l + c)xn-3 + B N 



(A19) 



and follows from ( A18[ ). 

Consequently, {X^) maximizes (minimizes) the 
normalization ratio Xn+i/xn for a given P and S, which 
proves (ii) and (Hi). The inequality (iv) then follows by 
taking the indicated limit, S->oo. 



[I] R. Healey, Stud. Hist. Philos. Sci. B, 

http://dx.doi.Org/10.1016/j.shpsb.2011.05.001 (2012). 

[2] C.M.S. Collaboration, Phys. Rev. Lett. 105, 032001 
(2010). 

[3] M.W. Zwierlein, C.A. Stan, C.H. Schunck, S.M.F. Rau- 
pach, S. Gupta, Z. Hadzibabic, and W. Ketterle, Phys. 
Rev. Lett. 91, 250401 (2003). 

[4] T. Bourdel, L. Khaykovich, J. Cubizolles, J. Zhang, 
F. Chevy, M. Teichmann, L. Tarruell, S.J.J.M.F. Kokkel- 
mans, and C. Salomon, Phys. Rev. Lett. 93, 050401 
(2004). 

[5] M. Combescot, O. Betbeder-Matibet, and F. Dubin, 

Phys. Rev. A 76, 033601 (2007). 
[6] C. K. Law, Phys. Rev. A 71, 034306 (2005). 
[7] C. Chudzicki, O. Oke, and W. K. Wootters, Phys. Rev. 



Lett. 104, 070402 (2010). 
[8] J. Schliemann, J.I. Cirac, M. Kus, M. Lewenstein, and 

D. Loss, Phys. Rev. A 64, 022303 (2001). 
[9] S. S. Avancini, J. R. Marinelli, and G. Krein, J. Phys. 

A: Math. Theor. 36, 9045 (2003). 
[10] S. Rombouts, D. V. Neck, K. Peirs, and L. Pollet, Mod. 

Phys. Lett. A17, 1899 (2002). 
[11] M. Combescot, O. Betbeder-Matibet, and F. Dubin, 

Phys. Rep. 463, 215 (2008). 
[12] M. Combescot and O. Betbeder-Matibet, Phys. Rev. 

Lett. 104, 206404 (2010). 
[13] M. Combescot, S.-Y. Shiau, and Y.-C. Chang, Phys. Rev. 

Lett. 106, 206403 (2011). 
[14] P. Sancho, J. Phys. A: Math. Theor. 39, 12525 (2006). 
[15] R. Ramanathan, P. Kurzynski, T.K. Chuan, M.F. San- 



9 



tos, and D. Kaszlikowski, Phys. Rev. A 84, 034304 

(2011) . 

[16] A. Gavrilik and Y. Mishchenko, Phys. Lett. A 376, 1596 

(2012) . 

[17] M. C. Tichy, F. Mintert, and A. Buchleitner, J. Phys. B: 
At. Mol. Opt. Phys. 44, 192001 (2011). 

[18] M. Combescot, X. Leyronas, and C. Tanguy, Europ. 
Phys. J. B 31, 17 (2003). 

[19] M. Combescot, Europhys. Lett. 96, 60002 (2011). 

[20] PA. Macmahon, Combinatory Analysis (Cambridge Uni- 
versity Press, Cambridge, 1915). 

[21] G. H. Hardy, J. E. Littlewood, and G. Polya, Inequalities 
(Cambridge University Press, Cambridge, 1988). 

[22] P. Harremoes and F. Tops0e, IEEE Trans. Inf. Theo. 47, 
2944 (2001). 

[23] K. Zyczkowski, Open Sys. Inf. Dyn. 10, 291 (2003). 



[24] J. Klotz, Techn. Rep. (University of Wisconsin) 591 
(1979). 

[25] T.-C. Wei, K. Nemoto, P.M. Goldbart, P.G. Kwiat, 
W.J. Munro, and F. Verstraete, Phys. Rev. A 67, 022110 
(2003). 

[26] A. G. Munford, Amer. Statist. 31, 119 (1977). 

[27] O. Giraud, J. Phys. A: Math. Theor. 40, F1053 (2007). 

[28] K. Zyczkowski, K. A. Penson, I. Nechita, and B. Collins, 
J. Phys. A: Math. Theor. 52, 062201 (2011). 

[29] D. N. Page, Phys. Rev. Lett. 71, 1291 (1993). 

[30] P. Hayden, D. W. Leung, and A. Winter, Comm. Math. 
Phys. 265, 95 (2006). 

[31] R. Ramanathan, A. Soeda, T. K. Chuan, and D. Kasz- 
likowski, arXiv:1108.2998vl (2011). 

[32] S. Mase, Ann. Inst. Statist. Math. 44, 479 (1992). 



